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The structure of the gluon propagator* 
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The gluon propagator has been calculated for quenched QCD in the Landau gauge at /3 = 6.0 for volumes 
16 3 x 48 and 32 3 x 64, and at f3 — 6.2 for volume 24 3 x 48. The large volume and different lattice spacings allow 
us to identify and minimise finite volume and finite lattice spacing artefacts. We also study the tensor structure 
of the gluon propagator, confirming that it obeys the lattice Landau gauge condition. 



1. Introduction 

The infrared behaviour of the gluon propaga- 
tor is important for an understanding of confine- 
ment. Previous conjectures range from a strong 
divergence to a propagator that vanishes in 
the infrared . Lattice QCD should be able 
to resolve this issue by first-principles, model- 
independent calculations. However, previous lat- 
tice studies [||J|] have been inconclusive since they 
have not been able to access sufficiently low mo- 
menta. Here we will report results using an asym- 
metric lattice with a spatial length of 3.3 fm This 
gives us access to momenta as small as 400 MeV. 

2. Lattice formalism 

We use the 'symmetric' definition of the gluon 
field, given by U^{x) = exp(igoaA^(x + A/2))- 
This gives the gluon field in momentum space, 



,-ig^a/2 

2ig$a 



B M (g)-~TrB M (<7) 



(1) 



where q denotes the discrete momenta = 
27ro M /(aL M ), B^q) ee U^q) - Ul(-q), and 
Ufj,(q) ee Yl x e ~ l ^ X Un( x )- An alternative, 'asym- 
metric' definition of the gluon field can be pro- 
vided by U/j,(x) — exp(igoaA' (x)). In momen- 
tum space, this differs from A fi (x) by a factor 
exp(ig (tl a/2). 

"Talk presented by J. I. Skullerud 
tUKQCD Collaboration 



Table 1 

Simulation parameters. The lattice spacing is 
taken from the string tension ||. 



Name 


a 


a- 1 (GeV) 


Volume 


Aconf 


Small 


6.0 


1.885 


16 a x 48 


125 


Large 


6.0 


1.885 


32 3 x 64 


75 


Fine 


6.2 


2.63 


24 3 x 48 


223 



The gluon propagator (q) is defined as 
Df u (q) = {A;(q)AU-q))/V, (2) 

where A^q) = t a A"(q). In the continuum Lan- 
dau gauge, the propagator has the structure 

DfM = 5 a \5^- qj ^f)D{q 2 ), 



(3) 



At tree level, D(q 2 ) will have the form D^°\q 2 ) = 
1/g 2 . On the lattice, this becomes D^\q) = 
a 2 / (4^2^s'm 2 (q ll a/2)). Since QCD is asymptoti- 
cally free, we expect that up to logarithmic cor- 
rections, q 2 D(q 2 ) — > 1 in the ultraviolet. Hence 
we define the new momentum variable q by q^ = 
(2/a) sm(q fJ- a/2), and use this throughout 

We have analysed three lattices, with different 
values for the volume and lattice spacing. The 
details are given in table |l|. All the configura- 
tions have been fixed to Landau gauge with an 
accuracy ((c^) 2 ) < lO" 12 . 
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3. Tensor structure 

By studying the tensor structure of the gluon 
propagator, we may be able to determine how 
well the Landau gauge condition is satisfied, and 
also discover violations of continuum rotational 
invariance. 

The continuum tensor structure (||) follows 
from the condition q^A^ = 0. This translates 
directly to the lattice provided we use the sym- 
metric definition of the gluon field. If we use the 
asymmetric definition, we will instead obtain the 
condition ^^sing^ + cos 9m — = 0- 

The tensor structure may be measured directly 
by taking the ratios of different components of 
Dfj_ u (q) for the same value of q. The results for 
the small lattice are summarised in table 0, and 
compared to what one would expect from (|3|) , and 
to what one would obtain by replacing q with q 
m (|). The results are similar for the two other 
lattices. It is clear from table || that our numerical 
data are consistent with the expectation from (||) . 
We can also see that in general, the asymmetric 
definition A' of the gluon field gives results which 
are inconsistent with this form. 

4. Finite size effects and anisotropics 

In the following, we are particularly interested 
in the deviation of the gluon propagator from the 
tree level form. We will therefore factor out the 
tree level behaviour and plot q 2 D(q 2 ) rather than 
D(q 2 ) itself. 

Fig. |l| shows the gluon propagator on the small 
lattice as a function of qa. For low momentum 
values on the small lattice, there are large dis- 
crepancies due to finite volume effects between 
points representing momenta along the time axis 
and those representing momenta along the spa- 
tial axes. These discrepancies are absent from 
the data from the large lattice, shown in fig. 
This indicates that finite volume effects here are 
under control. 

However, at higher momenta, there are 
anisotropics which remain for the large lat- 
tice data, and which are of approximately the 
same magnitude for the two lattices. These 
anisotropics are considerably reduced on the fine 
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Figure 1. The gluon propagator multiplied by q 2 
as a function of q for the small lattice. The filled 
triangles denote momenta directed along the time 
axis, while the filled squares denote momenta di- 
rected along one of the spatial axes. 
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Figure 2. The gluon propagator multiplied by 
q 2 as a function of q for the large lattice. The 
symbols are as in fig. [I]. 

lattice, indicating that they arise from finite lat- 
tice spacing errors. In order to eliminate these 
anisotropies, we select momenta lying within a 
cylinder of radius Aqa = 2 x 27r/32 along the 4- 
dimcnsional diagonals. The result of this cut on 
the large lattice is shown in fig. ||. A more de- 
tailed discussion of these cuts can be found in . 

5. Conclusions 

We have evaluated the gluon propagator on 
three different lattices. The tensor structure has 
been analysed and shown to agree with the con- 
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Table 2 

Tensor structure for the small lattice, q is in units of 2ir/L s , where L s is the spatial length of the lattice. 
The theoretical predictions are the values for the ratios one obtains from (||), and from (||) with q — > q. 
The numbers in brackets are the statistical uncertainties in the last digit(s). Where no error is quoted, 
the statistical uncertainty is less than 10 -6 . 



Theoretical prediction This simulation 



Q 


Components 


Using q 


Using q 


Using A 


Using A' 


[2,1,0,0] 


(M)/(l,2) 


-0.5 


-0.509796 


-0.509796 


-0.519783 




(M)/(2,2) 


0.25 


0.259892 


0.259892 


0.259892 




(l,l)/(3,3) 


0.2 


0.206281 


0.204(8) 


0.204(8) 




(l,2)/(2,2) 


-0.5 


-0.509796 


-0.509796 


-0.5 


[4,1,0,0] 


(1,1)/(1,2) 


-0.25 


-0.275899 


-0.275899 


-0.331821 




(l,l)/(3,3) 


0.05882 


0.070736 


0.076(3) 


0.076(3) 




(l,2)/(2,2) 


-0.25 


-0.275899 


-0.275899 


-0.229402 


[4,2,1,0] 


(1,1)/(1,2) 


-0.625 


-0.681848 


-0.678(9) 


-0.743(10) 




(l,l)/(2,3) 


-2.5 


-2.47137 


-2.3(4) 


-2.5(5) 




(l,2)/(2,2) 


-0.4706 


-0.502914 


-0.500(6) 


-0.456(6) 
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Figure 3. The gluon propagator multiplied by q 2 
as a function of q for the large lattice, after the 
cylindrical cut. 



tinuum Landau gauge form. By studying the 
anisotropics in the data, we have been able to 
conclude that finite volume effects are under con- 
trol on the largest of our lattices. 

A clear turnover in the behaviour of q 2 D(q 2 ) 
has been observed at q ~ lGeV, indicating that 
the gluon propagator diverges less rapidly than 
1/q 2 in the infrared, and may be infrared fi- 
nite. An analysis of scaling and the functional 
behaviour of these results is presented in M . 
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